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FOREWORD 

This  report  was  prepared  under  Work  Order  No.  086412  of  Contract  No. 
N00014-78-C-0204  in  support  of  a,  U.S.  Naval  Postgraduate  School  research 
project  sponsored  by  the  Naval  Air  Systems  Command,  AIR  370,  and  Naval 
Avionics  Center,  Indianapolis.  This  report  describes  the  results  of  the 
second  stage  analysis  of  aircraft  measurements  of  micrometeorological  param¬ 
eters  made  at  Panama  City,  Florida,  in  1978.  ^Thg_j£trs±_stage  of  analysis 
has  been  given  in  a  previous  report  (Technical  Report  No.  BDM/M- 008-79}^ 

The  work  was  done  in  cooperation  with  Dr.  Ralph  Markson  and  Mr.  Jan  Sedlacek 
of  Airborne  Research  Associates  and  Drs.  Ken  Davidson  and  Gordon  Schacher  of 
the  Naval  Postgraduate  School. 
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ABSTRACT 

-"''The  dissipation  rate  of  turbulent  kinetic  energy,  e,,  and  the  tempera- 
ture  structure  function  parameter,  vCj  ,  -have-  been  measured  over  water  from 
the  near  surface  (Z  =  3m)  to  the  top\of  the  boundary  layer.  The  near  sur- 
face  values  of  and  C^were  used  to  calculate  the  velocity  and  temperature 

Monin-Obukhov  scaling  parameter$,.u*_afld--T^>  The  unstable  data  were  used  to 

t  2 ' 

evaluate  the  feasibility  of  extrapolating  the  values  of  js)and  C-j-t  as  a  func¬ 
tion  of  height  with  empirical  scaling  formulae.  The  dissipation^ rate  scaling 

formula  of  Wynqaard  et  al.  (1971a)  qave  a  very  qood  fit /to  an  averaqe  of  the 
.y*  <•?  ?  >  ■  J  -  -4/3  '- 

e  data  for  Z  <  I.  -  100.  Although  the  /data  obeyed  a [1  J'  height  depen-  - 

dence  for  0.1  £>,<  Z  <  0.5  Z>>»  the  scalinq  formula  of  Wynqaard  et  al.  (1971b) 

disagreed  with  the  average  Cy  by  as  much  as  50%.  At  this  point  it  is  not 

clear  if  the  discrepancy  is  a  unique  property  of  the  marine  boundary  layer 

or  if  it  is  simply  some  unknown  instrumental  or  analytic  problem. 


THE  BDM  CORPORATION 


A.  INTRODUCTION 


This  paper  is  a  report  on  measurements  of  temperature  structure  param- 
2 

eter,  Cy  ,  and  rate  of  dissipation  of  turbulent  kinetic  energy,  e,  made  from 
a  light  aircraft  using  microthermal  sensors  as  part  of  a  study  of  turbulence 
in  the  marine  boundary  layer.  A  typical  flight  included  measurements  as  low 
as  three  meters  above  the  sea  surface  in  an  effort  to  characterize  the  sur¬ 
face  layer  turbulence  parameters.  The  flights  were  usually  extended  well 
above  the  inversion  -  a  maximum  altitude  of  three  kilometers  being  typical. 

The  data  were  gathered  near  Panama  City,  Florida,  as  part  of  a  program 
of  marine  boundary  layer  research,  the  ultimate  goal  of  which  is  the  formu¬ 
lation  of  a  model  that  will  allow  reasonable  estimates  of  atmospheric  turbu¬ 
lence  and  mean  properties  from  the  surface  throughout  the  entire  marine 
boundary  layer  and  immediately  above.  Such  a  model  would  have  application 
to  such  varied  subjects  as  turbulent  transport  of  heat,  water  vapor  and 
pollutants,  optical  propagation,  aerosols,  radar  propagation,  and  communi¬ 
cations. 

2 

For  example,  Cj  is  important  for  optical  propagation  studies  due  to 

2 

its  relation  to  the  index  of  refraction  structure  parameter,  CN  ,  (Friehe, 
1977), 

CN2  =  (79  x  10'6  P/T2)2  (Cy2  +  .11  CTQ  +  3.2  x  10-3  CQ2)  (1) 

2 

where  Cp  is  the  humidity  structure  parameter,  Cyp  the  cospectrum  structure 
parameter,  P  is  the  pressure  in  mb,  T  the  absolute  temperature,  and  Q  is  in 
gm/m  .  Overland  the  water  vapor  fluctuation  contribution  is  usually  negli¬ 
gible,  so  one  can  write 

CN2  =  (79  x  10"6  P/T2)2  Cy2  (2) 

2  2 

This  relationship  is  shown  in  Figure  1  with  /Cy  as  a  function  of  altitude 
for  the  U.S.  Standard  Atmosphere.  In  the  marine  boundary  layer  equation  (2) 
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is  only  an  approximation  because  of  the  large  vertical  water  vapor  fluxes 

that  can  be  present  over  the  ocean  (Wesley,  1976).  Friehe  et  al.  (1975) 

2  2 

found  for  the  San  Diego  FLIP  data  that  the  Cy  contribution  to  CN  was 
typically  70%.  Ideally,  one's  model  should  be  equally  adept  at  handling 
temperature  and  water  vapor  height  dependences. 

This  paper  will  focus  on  the  utilization  of  the  data  to  examine  certain 
aspects  of  the  Monin-Obukhov  similarity  (MOS)  expressions  for  the  dimension¬ 
less  temperature  structure  function  and  the  dimensionless  dissipation  func¬ 
tion  in  the  marine  atmospheric  surface  layer.  The  validity  of  the  existing 
empirical  expressions  of  the  asymptotic  height  dependences,  variance  of  the 
data  from  surface  extrapolations,  and  the  influence  of  the  height  of  the 
boundary  layer  will  all  be  explored.  Of  the  nine  days  on  which  flights  were 
made,  seven  were  characterized  by  unstable  surface  layers  and  two  by  stable 
surface  layers.  Since  the  height  of  the  boundary  layer  was  so  low  on  the 
stable  days  (typically  100  meters),  the  influence  of  the  height  of  the 
boundary  layer  tended  to  dominate  the  surface  layer  conditions  causing  most 
of  the  quantitative  analyses  mentioned  immediately  above  to  be  of  dubious 
relevance.  In  view  of  this,  the  analysis  of  the  surface  layer  scaling 

properties  was  restricted  to  the  unstable  data. 

B.  THEORETICAL  BACKGROUND 

The  boundary  layer  is  that  part  of  the  atmosphere  where  friction  with, 
and  heating  by  the  surface  play  an  important  part  in  the  generation  of  tur¬ 
bulence.  Near  the  surface  the  turbulence  properties  can  be  scaled  in  terms 

of  the  Reynolds  stress,  tq,  and  the  surface  sensible  virtual  temperature 
flux,  Q0y,  using  Monin-Obukhov  similarity,  where  departures  from  neutral 
equilibrium  are  parameterized  by  Z/L.  Those  interested  in  a  more  complete 
treatment  of  the  surface  layer  should  refer  to  Businger  (1973)  and  Wyngaard 
et  al .  (1971a  and  1971b).  The  velocity  and  temperature  are  scaled  by  u*  and 
T*  and  stability,  L, 
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t0  =  p  u* 


"U*Ty* 


L  "  ^  W 

where  p  is  the  density  of  air,  T  is  the  absolute  temperature,  Ty  the  virtual 
potential  temperature,  g  the  acceleration  of  gravity  and  k  is  von  Karmon's 
constant  (we  have  used  k  =  0.35).  Ty*  is  related  to  the  water  vapor  mixing 
ratio  scaling  parameter,  q*,  by 

Ty*  =  T*  +  6.1  x  10'4  T  q*  (4) 

where  q  is  given  in  gm/kg. 

In  the  inertial  subrange  of  locally  isotropic  turbulence,  fluctuations 
in  horizontal  wind  speed  and  temperature  can  be  characterized  by  a  one¬ 
dimensional  fourier  power  spectrum  S(k)  (Tatarski,  1961) 

Su(k)  =  0.52  e2/3  k"5/3  (5a] 

ST(k)  =  0.25  CT2  k"5/3  (5b) 

where  k  is  the  wavenumber,  defined  such  that 

oo 

/  Sy(k)  dk  =  < ( x  -  x)2>  (6) 

o 

The  values  of  various  empirical  constants  (such  as  k)  were  taken  from 
Champagne  et  al .  (1977). 

In  the  surface  layer,  one  can  write  expressions  for  e  (Wyngaard  et  al . , 
1971a)  and  Cy2  (Wyngaard  et  al . ,  1971b) 
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e=^*e(5)  (7a) 

CT2  =  T*2  Z“2/3  f(c)  (7b) 

where  £  =  Z/L.  The  present  form  of  the  dimensionless  dissipation  function, 
<f>e  (5),  is 

*e  (C)  =  (1  +  0.5  U|2/3)3/2  5  <  0  (8a) 

*£  (C)  =  (1  +  2.5  \K\2/3)3/2  5  >  0  (8b) 

while  that  for  the  dimensionless  structure  functions  parameter,  f(s),  is 

f (5)  =  4.9  (1  -  7s)"2/3  5  <  0  (9a) 

f (C)  =  4.9  (1  +  2.4  K  2/3)  5  >  0  (9b) 

The  principle  of  similarity  is  invoked  under  the  assumption  that  |L|  is 


considerably  less  than  the  height  of  the  boundary  layer,  h.  As  was  previously 
mentioned,  this  condition  was  not  met  for  the  stable  data.  Under  unstable 
conditions  the  boundary  layer  height  is  determined  by  the  height  of  the  lowest 
inversion,  Z.  (Deardorff,  1974),  so  that  h  =  Z^. 

C.  INSTRUMENTATION  AND  ANALYSIS 


The  platform  for  these  measurements  is  a  single  engine  turbo-charged 
Bellanca  Viking  operated  by  Airborne  Research  Associates.  The  aircraft  is 
well  instrumented  allowing  simultaneous  measurements  of  air  temperature, 
pressure,  dew  point,  electric  field,  and  infrared  surface  temperature. 
Details  of  the  instrumentation  and  data  acquisition  are  contained  in  a  tech¬ 
nical  report  (Fa i rail  and  Markson,  1979).  We  will  only  discuss  several 
aspects  of  the  turbulence  measurements  deemed  relevant  to  the  interpretation 
of  the  data. 
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Dissipation  rate,  e,  was  determined  from  velocity  fluctuations  sensed 
with  a  constant  temperature  anemometer  employing  tungsten  wires  4.5  p  in 
diameter  and  1.7  mm  in  length  operated  at  a  50%  overheat.  The  signal  was 
bandpass  filtered  (f  =  50  Hz  and  f  =  500  Hz)  before  being  squared  and 

X  U 

averaged.  One  can  calculate  e  from  the  filtered  mean  square  fluctuations 
using 


where  u  =  60  m/sec  is  the  average  aircraft  true  airspeed. 

2 

Temperature  structure  function,  Cy  ,  was  determined  from  temperature 
difference  fluctuations  sensed  with  pairea  tungsten  wires  separated  a  dis¬ 
tance  d  =  0.85  meters.  Two  complete  systems  were  used  (one  ac  Wheatstone 
bridge  and  one  dc)  each  employing  an  independent  pair  of  wires.  The  two 
systems  systematically  disagreed  by  about  20%,  possibly  due  to  slightly 
different  frequency  response  characteristics.  The  velocity  sensitivity  of 
the  sensors  as  a  function  of  bridge  current  was  carefully  measured  at 
IF  =  60  m/sec  in  the  laboratory  to  optimize  the  temperature  fluctuation 
measurements.  The  mean  square  temperature  fluctuation  data  were  corrected 
for  velocity  effects  using  the  known  velocity  sensitivity  and  the  values  of 
velocity  fluctuation  intensity  measured  simultaneously  with  the  hot-wire 
anemometer.  Furthermore,  one  temperature  bridge  was  operated  at  half  the 
sensor  current  of  the  other  in  order  to  flag  possible  velocity  sensitivity 

problems,  should  they  occur.  In  practice,  the  velocity  contribution  to  the 

2 

temperature  signal  was  considerably  less  than  the  system  noise.  Cy  was 
calculated  using  the  equation 


=  <{T  -  T)2> 


d 


d-2/3 


(ID 


which  is  based  on  the  assumption  that  the  probe  separation,  d,  is  within  the 
inertial  subrange. 
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Because  the  larger  scale  eddies  become  more  restricted  as  one  approaches 

2 

the  sea  surface,  the  measured  values  of  Cy  and  e  are  subject  to  error  due 
to  the  assumptions  of  isotropy  in  eqs.  10  and  11.  This  is  particularly 
important  for  this  data  because  of  the  unusually  low  altitudes  (Z  =  3  meters) 
employed.  We  have  examined  this  question  by  replacing  the  inertial  subrange 
forms  of  Su  and  Sy  with  forms  valid  over  the  entire  range  of  k.  Thus,  we 
can  compute  the  actual  measured  values  of  the  mean  square  fluctuations 


<(u  -  u>\u  *  / 


Su(f)  d  f 


(12a) 


<(T  -  T)2>d  =2  J  [1  -  cos  (kd)]  Sy(k)  d  k  (12b) 

o 

where  we  have  used  the  formulae  for  Su  and  Sy  given  by  Kaimal  et  al . ,  (1972) 

for  neutral  stability.  Since  Kaimal  et  al . ,  (1972)  indicate  that  the  size 

scales  for  unstable  conditions  are  considerably  larger  than  the  neutral  case, 

this  calculation  should  serve  to  establish  an  upper  limit  on  the  effects. 

We  have  expressed  the  result  in  terms  of  the  ratio  of  the  measured  to  the 

2  2 

actual  value  (em/e  or  Cym  /Cy  )  as  a  function  of  altitude  (Figure  2)  for 

those  filter  frequencies  and  probe  spacing  used  for  the  experiment.  Based 

on  this  calculation,  one  can  conclude  that  the  surface  effects  on  the  measure- 

2 

ments  of  e  were  negligible  and  that  the  upper  limit  of  the  effects  on  the  Cy 
were  on  the  order  of  10%. 


D.  FLIGHT  PATTERN 


Aircraft  are  particularly  useful  for  atmospheric  research  because  they 
are  able  to  sample  large  pieces  of  atmosphere  in  relatively  short  time  periods. 
Unfortunately,  even  at  sixty  meters  per  second,  one  cannot  expect  to  transverse 
the  larger  scale  eddies  normally  associated  with  turbulence  in  just  a  few 
seconds.  As  a  result,  if  one  wishes  to  define  the  turbulence  properties  of 
the  atmosphere  to  some  reasonable  statistical  confidence,  an  averaging  period 
of  a  few  minutes  (on  the  order  of  10  kilometers)  is  required.  Obviously  for 
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Figure  2.  The  Expected  Effect  of  the  Lower  Isotropic  Limit  on  the  Turbulence 
Measurements  as  the  Aircraft  Nears  the  Sea  Surface 

The  measurement  error  is  expressed  as  the  ratio  of  the  measured  value  to  the 
actual  value  ( Xm/X)  for  (dashed  line)  and  Cj2  (solid  line).  The  short  hori¬ 
zontal  line  with  the  vertical  arrow  indicates  the  minimum  altitude  of  the  air¬ 
craft. 
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|  the  surface  layer  investigation  one  not  only  requires  very  low  altitude  data 

but  a  greater  density  of  measurement  altitudes  near  the  surface.  The  final 
|  flight  pattern  decided  upon  consists  of  a  series  of  constant  altitude  runs 

of  two  minutes  duration  each  at  successively  greater  altitude.  In  order  to 
I  maintain  the  same  approximate  location,  the  flight  direction  was  reversed 

*  at  alternate  altitudes.  Normally  the  process  was  started  at  the  lowest 

r  altitude  of  three  meters,  as  "eyeball"  estimated  by  the  pilot.  At  this  point 

.  the  aircraft  altimeter  was  set  to  ten  feet  and  subsequent  altitudes  were 

based  on  this  setting.  The  altitudes  of  the  level  runs  were  increased  approx¬ 
imately  exponentially. 

This  type  of  flight  pattern  has  been  named  a  "ladder  profile"  with  the 
two  minute  runs  being  the  rungs  of  the  ladder.  The  total  time  required  to 
execute  a  complete  ladder  profile  was  on  the  order  of  30  minutes. 


i 

I 

I 

I 

I 

l 

I 

I 

I 

I 
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E.  SURFACE  LAYER  SCALING  PARAMETERS 


The  Monin-Obukhov  scaling  parameters  for  each  profile  were  obtained  from 
the  lowest  two  or  three  data  points.  First,  q*  was  calculated  from  a  simple 
bulk  formulation  using  a  ten  meter  drag  coefficient 

q*  =  0.036  (q  -  q$)  (13) 

where  q  is  the  mixing  ratio  at  2  =  10  meters  and  qs  is  the  value  for  the  sea 

surface  calculated  from  the  sea  surface  temperature  assuming  100%  relative 

humidity.  Given  the  relatively  smaller  contribution  of  q*  to  L,  a  more 

sophisticated  treatment  of  water  vapor  was  not  considered  to  be  worthwhile. 

At  this  point,  values  of  u*  and  T*  were  selected  that  gave  good  fits  of 

eqs.  7a  and  7b  to  the  lower  altitude' values  of  e  and  Cy  .  Since  |c|  was 

usually  much  less  than  unity  in  this  region,  the  actual  form  of  6  (t)  and 

£ 

f(5)  would  be  de-emphasized  in  the  determination  of  u*,  T*,  and  L.  The  sign 
of  T*  was  determined  from  the  air-sea  temperature  difference.  The  scaling 
parameters  for  all  twenty-one  profiles  taken  during  PC  II  are  given  in  Table  1. 
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TABLE  1.  SURFACE  LAYER  SCALING  PARAMETERS  FOR  THE  LADDER  PROFILES 
|  OBTAINED  DURING  PC  II 

N  is  the  number  of  points  in  each  profile;  the  other  quantities  are  defined 
I  in  the  text. 


Profile 

Date 

Time 

u* 

T* 

L 

h 

N 

# 

(m/s) 

(°C) 

(gm/kg) 

(m) 

(m) 

1 

11/26 

1252 

0.400 

-O.OS2 

-0.160 

-125 

850 

14 

2 

11/26 

1436 

0.230 

-0.095 

-0.160 

-  50 

900 

14 

3 

12/02 

1405 

0.240 

-0.135 

-0.180 

-  29 

230 

18 

4 

12/03 

1108 

0.290 

0.030 

0.000 

233 

200 

11 

5 

12/03 

1201 

0.240 

0.015 

-0.018 

420 

60 

10 

6 

12/03 

1232 

0.180 

0.024 

-0.009 

120 

60 

9 

7 

12/03 

1339 

0.340 

0.031 

-0.040 

403 

120 

11 

8 

12/05 

1532 

0.240 

-0.260 

-0.410 

-  15 

700 

10 

9 

12/05 

1624 

0.260 

-0.270 

-0.400 

-  16 

700 

13 

10 

12/07 

1511 

0.250 

0.015 

-0.007 

390 

75 

12 

11 

12/07 

1601 

0.235 

0.025 

0.010 

171 

100 

11 

12 

12/10 

1259 

0.530 

-0.440 

0.000 

-  53 

760 

10 

13 

12/10 

1324 

0.380 

-0.350 

0.000 

-  34 

760 

11 

14 

12/10 

1410 

0.320 

-0.490 

-0.490 

-  15 

760 

15 

15 

12/10 

1523 

0.340 

-0.480 

-0.470 

-  17 

820 

14 

16 

12/10 

1637 

0.340 

-0.490 

-0.500 

-  16 

980 

14 

17 

12/11 

1021 

0.280 

-0.440 

-0.430 

-  13 

700 

12 

18 

12/12 

1642 

0.280 

-0.180 

-0.500 

-  24 

800 

18 

19 

12/13 

1154 

0.190 

■  -0.210 

-0.470 

-  10 

600 

17 

20 

12.13 

1459 

0.170 

-0.200 

-0.420 

-  9 

550 

14 

21 

12.13 

1721 

0.140 

-0.120 

-0.440 

-  5 

450 

17 

1 

i 

I 
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F.  LADDER  PROFILES 


The  locations  of  the  profiles  are  indicated  in  Figure  3  by  profile  num¬ 
ber.  Profiles  12  and  13,  which  were  overland,  were  not  used  in  the  results 
presented  in  Section  7.  Since  the  profile  data  are  completely  cataloged 
elsewhere  (Fai rail  and  Markson,  1979),  we  will  give  only  a  few  examples  here. 

The  situation  encountered  for  the  stable  surface  layers  is  illustrated  in 

2 

Figure  4  where  the  deviations  of  Cy  from  the  surface  layer  expression 
(eq.  7b)  extend  down  to  the  lowest  altitudes.  The  example  given  is  one  of 
the  extremes.  There  were  other  stable  profiles  where  the  fit  was  much  better, 
but,  in  general,  the  stable  surface  layers  appeared  to  be  dominated  by  the 
shallowness  of  the  boundary  layer. 

The  boundary  layer  was  much  more  extensive  on  unstable  days  (the  inver¬ 
sion  height  was  typically  800  meters).  In  profile  1  (Figure  5a)  the  surface 

2 

layer  expression  for  Cy  is  a  good  fit  up  to  at  least  300  meters.  Note  that 
in  this  case  the  upper  altitude  values  of  e  are  smaller  than  expected  based 
on  the  near  surface  values,  while  in  profile  19  (Figure  5b)  the  reverse  is 
true.  This  effect  may  be  due  to  the  fact  that  we  ignored  the  possible  exis¬ 
tence  of  wind  roll  vortices  (LeMone,  1973)  or  other  two-dimensional  structure. 
2 

In  Figure  5b  the  CT  data  displays  a  feature  that  often  occurred.  The  higher 

-4/3 

altitude  values  are  very  well  fit  by  a  Z  '  curve,  but  the  actual  values  are 

2 

not  consistent  with  the  near  surface  values,  the  near  surface  values  of  CT 

-4/3 

being  somewhat  lower  than  expected  on  the  basis  of  the  Z  region.  This 
aspect  will  be  discussed  in  Section  7  and  Section  8. 

G.  SURFACE  LAYER  DIMENSIONLESS  TURBULENCE  RESULTS 

Given  values  of  u*  and  T*  for  an  unstable  profile,  the  dimensionless 
3  2  2/3  2 

quantities  e  k  Z/u*  and  Cy  Z  '  /(4.9  T*  )  were  calculated  for  each  data 

point.  Because  u*  and  T*  were  determined  from  the  near  surface  values  of  e 
2 

and  Cy  ,  the  dimensionless  quantities  automatically  have  values  near  unity 
for  -Z/L  «  1.  The  dimensionless  values  were  then  averaged  in  bins  of  diff¬ 
erent  dimensionless  stability  length  allocated  according  to  log|Z/L|.  Since 


The  data  points  plotted  are  virtual  potential  temperature  (+),  dew  point  temperature  (*),  Cj?  (*)» 
and  e(o).  The  solid  line  is  the  MOS  expression  for  Cj2  and  the  long-dash  line  is  the  MOS 
expression  fore.  The  extreme  left-hand  side  of  the  graph  shows  an  expanded  plot  of  virtual  poten 
tial  temperature. 
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the  number  of  profiles  taken  on  a  given  flight  varied  from  one  to  five,  a 
weighted  average  was  used  to  avoid  overemphasizing  the  conditions  that  happen 
to  exist  during  those  flights  with  more  profiles.  Each  profile  was  assigned 
a  weight  W  =  1/M  where  M  was  the  number  of  overwater  profiles  in  that  flight. 
In  order  to  eliminate  effects  due  to  the  inversion,  the  data  had  to  be  res¬ 
tricted  based  upon  .  The  criteria  that  allowed  the  most  data  to  be  used 
without  introducing  inversion  effects  were 


Z  <  Z.  -  100  for  e 

(14a) 

Z  <  0.5  Z.  for  CT2 

(14b) 

The  results  are  given  in  Figures  6a  and  6b  with  the  solid  curves  representing 
eq.  8a  and  eq.  9a,  respectively. 

In  order  to  examine  the  absolute  height  dependence  (rather  than  the 
dimensionless  height  dependence)  the  quantity  called  R  was  calculated  for 
each  data  point.  R  is  the  ratio  of  the  dimensionless  turbulence  value  to 
that  value  one  would  expect  based  on  eq.  8a  and  9a. 

R  =  (1  +  0.5|Z/L|2/3)‘3/2  (15a) 

E  u* 

C  2  Z2/3 

Rt  =  — - r  (1  *  7  Z/L)2/3  (15b) 

4.9  V 

The  values  of  R  were  then  averaged  in  bins  of  different  altitude  in  a  manner 
similar  to  that  just  described.  The  averages  of  R  are  shown  in  Figure  7,  and 
the  ratio  of  the  standard  deviation,  aD,  to  R  is  shown  in  Figure  8.  Taken 
together,  these  two  figures  illustrate  one's  ability  to  predict  e  and  Cy 
as  a  function  of  height  from  near  surface  measurements.  The  20 %  standard 
deviation  for  the  near  surface  values  (Figure  8)  is  due  to  the  scatter 
inherent  in  the  two  minute  averages. 
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100 


-Z/L 


Figure  6a.  Dimensionless  Dissipation  Rate  as  a  Function  of 
Dimensionless  Atmospheric  Stability 

The  points  represent  the  weighted  average  of  the  profile  data,  the  vertical 
bars  represent  the  statistical  error  in  the  mean  estimate,  and  the  number  is 
the  weighted  number  of  points  in  the  average.  The  solid  line  is  from  eq.  8a. 
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Figure  6b.  Dimensionless  Temperature  Structure  Function  of 
Dimensionless  Atmospheric  Stability 


The  points  represent  the  weighted  average  of  the  profile  data,  the  vertical 
bars  represent  the  statistical  error  in  the  mean  estimate,  and  the  number  is 
the  weighted  number  of  points  in  the  average.  The  solid  line  is  from  eq.  9a 
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H.  DISCUSSION 


Kaimal  et  al . ,  (1976)  divide  the  unstable  boundary  layer  into  three 
regions:  the  surface  layer  (Monin-Obukhov  scaling),  the  free  convection 
layer,  and  the  mixed  layer.  Panofsky  (1977)  points  out  that  the  free  convec¬ 
tion  layer  is  a  region  where  Monin-Obukhov  scaling  and  mixed  layer  scaling 
both  apply.  Hence,  he  has  called  this  region  the  convective  matching  layer 
defined  by  -L  <  Z  <  0.1  Z.  At  this  point,  we  wish  to  consider  how  well  the 
data  fits  the  surface  layer  equations  and  the  altitude  limits  of  the  fit. 

Certainly  the  average  fit  of  the  e  data  to  eq.  8a  (Figures  6a  and  7)  is 
surprisingly  good.  Furthermore,  the  formula  gives  reasonable  results  through¬ 
out  the  majority  of  the  boundary  layer  (Z  <  Z^  -  100).  The  free  convection 
limit  of  eq.  7a  yields 


e  -  A  (g/T)  QqV,  -  Z/L  »  1  (16) 

where  A  =  0.36.  Lenschow  (1974)  has  found  A  *  0.43,  while  Kaimal  et  al. , 
(1976)  measured  values  of  A  between  0.5  and  0.7  overland.  Actually,  this 
limit  was  reached  only  on  profile  21  which  happened  to  be  more  consistent 
with  A  =  0.5. 

.  .  .  2 

In  Figure  6b  we  can  see  considerable  deviations  of  the  CT  from  eq.  9a. 
The  fact  that  the  data  lie  on  the  expected  curve  for  -Z/L  «  1  is  not  signi¬ 
ficant  since  the  values  of  T*  were  selected  to  force  a  fit  in  this  region. 

It  is  of  interest  to  speculate  on  the  possible  significance  of  this  disagree¬ 
ment.  Recall  that  eq.  9a  was  based  on  overland  measurements  (Wyngaard  et  al . , 
1971b)  with  only  a  few  points  for  -Z/L  >  1.  However,  the  more  recent  measure¬ 
ments  of  Kaimal  et  al . ,  (1976)  were  examined  in  the  convective  limit  of  eq.  7a 
(note  Qq  is  the  sensible  temperature  flux,  -u*  T*) 

CT2  =  2.67  (g/T)"2/3  QoV"2/3  Qq2  Z‘4/3  (17) 

and  excellent  agreement  was  found  for  0.1  Z^  <  Z  <  0.5  Z..  Thus,  one  can 
have  confidence  in  eq.  9a  for  the  overland  case.  The  overwater  measurements 
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of  Davidson  et  al . ,  (1979)  show  disagreement  with  eq.  9a  similar  to  that  of 
Figure  6b,  but  this  disagreement  is  attributed  to  probable  sea-salt  contamin¬ 
ation  of  the  sensors.  Frisch  and  Ochs  (1975)  found  deviations  in  the  Z“4^3 

height  dependence  of  eq.  17  overwater  for  Z  >  0.1  Z. .  However,  they  deter- 

-2/3  2  2  ' 

mined  QQy  '  Qq  from  values  of  Cy  measured  at  Z  =  30  m  and  did  not  look 

at  lower  altitudes.  A  corresponding  analysis  of  the  PC  II  data  (determining 
Q0v  Q0  from  the  value  of  Cy  at  Z  =  30  m)  shows  no  deviation  from  eq.  17 
for  0.1  Z.j  <  Z  <  0.5  Z..  The  greater  height  dependence  of  R  at  lower  alti¬ 
tudes  (Figure  7)  suggests  that  the  discrepancy  is  some  type  of  surface  effect. 

2 

Recall  that  in  Section  3  the  surface  effect  on  the  Cy  measurement  process 
was  estimated  to  have  an  upper  limit  of  10«  at  Z  =  3  m  (Figure  2).  Certainly 
this  does  not  preclude  the  possibility  of  some  other  measurement  error  such 
as  frequency  response  losses  or  saturation  effects,  but  at  this  writing  no 
such  effects  have  been  discovered.  Since  we  are  unable  to  propose  a  physical 
mechanism  for  this  discrepancy  in  the  marine  boundary  layer,  we  prefer  to 
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